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1 Introduction 

The Einstein Equivalence Principle (EEP) is a founding principle of relativ¬ 
ity [1]. One of the constituent elements of EEP is Local Lorentz Invariance 
(LLI), which postulates that the outcome of a local experiment is indepen¬ 
dent of the velocity and orientation of the apparatus. The central importance 
of this postulate has motivated tremendous work to experimentally test LLI. 
Also, a number of unification theories suggest a violation of LLI at some level. 
However, to test for violations it is necessary to have an alternative theory 
to allow interpretation of experiments [1], and many have been developed 
[2, 3, 4, 5, 6, 7]. The kinematical frameworks (RMS) [2, 3] postulate a sim¬ 
ple parameterization of the Lorentz transformations with experiments setting 
limits on the deviation of those parameters from their values in special relativ¬ 
ity (SR). Because of their simplicity they have been widely used to interpret 
many experiments [8, 9, 10, 11]. More recently, a general Lorentz violating ex¬ 
tension of the standard model of particle physics (SME) has been developed 
[6] whose Lagrangian includes all parameterized Lorentz violating terms that 
can be formed from known fields. 

This work analyses rotating laboratory Lorentz invariance experiments 
that compare precisely the resonant frequencies of two high-Q factor (or high 
finesse) cavity resonators. High stability electromagnetic oscillatory fields are 
generated by implementing state of the art frequency stabilization systems 
with the narrow line width of the resonators. Previous non-rotating experi¬ 
ments [12, 10, 13] relied on the rotation of the Earth to modulate putative 
Lorentz violating effects. This is not optimal for two reasons. Firstly, the sensi¬ 
tivity to Lorentz violations is proportional to the noise of the oscillators at the 
modulation frequency, typically best for periods between 10 and 100 seconds. 
Secondly, the sensitivity is proportional to the square root of the number of 
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periods of the modulation signal, therefore taking a relatively long time to 
acquire sufficient data. Thus, by rotating the experiment the data integra¬ 
tion rate is increased and the relevant signals are translated to the optimal 
operating regime [14]. 

In this work we outline the two most commonly used test theories (RMS 
and SME) for testing LLI of the photon. Then we develop the general frame¬ 
work of applying these test theories to resonator experiments with an emphasis 
on rotating experiments in the laboratory. We compare the inherent sensitivity 
factors of common experiments and propose some new configurations. Finally 
we apply the test theories to the rotating cryogenic experiment at the Univer¬ 
sity of Western Australia, which recently set new limits in both the RMS and 
SME frameworks [15]. Note added: Two other concurrent experiments have 
set some similar limits [16, 17]. 


2 Common Test Theories to Characterize Lorentz 
Invariance 

The most famous test of LLI (or the constancy of the speed of light) was that 
conducted by Michelson and Morley in 1887 [18] with a rotating table and a 
Michelson interferometer. In actual fact, the theoretical framework used by 
Michelson and Morley was not a test of LLI, since the concept did not exist 
at the time, but that of an aether drift. The relative motion of the apparatus 
through the aether was thought to induce a phase difference between the 
arms of the interferometer (and hence an interference pattern) depending 
on the orientation. Thus, as the Earth moved from one end of its orbit to 
the opposite end, the change in its velocity should be a detectable value. 
Michelson and Morley found no fringe shifts due to Earth motion around the 
sun and reported a null result. Since the Michelson Morley experiment, there 
have been many other types of experiments devised to test the validity of SR 
and the constancy of light. However, to interpret these experiments one must 
formulate an alternative test theory, and in this section we outline two of the 
most commonly used. 

2.1 Robertson, Mansouri, Sexl Framework 

A simple kinematic test theory that has been widely used is that of Robert¬ 
son, Mansouri and Sexl (RMS)[2, 3], where time standards (“clocks”) and 
length standards (“rods”) are considered without taking into account their 
underlying structure. This framework postulates a preferred frame S(T, X) 
which satisfies LLI, and a moving frame S(t, x), which does not. The prime 
candidate for the preferred frame is taken as the Cosmic Microwave Back¬ 
ground (CMB), since any anticipated non-symmetries are expected to arise 
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from Planck-scale effects during the creation of the universe. In this frame¬ 
work we analyse the Poynting vector direction of the electromagnetic signal 
with respect to the velocity of the lab through the CMB. 

The normal Lorentz Transformations for a boost in the x direction are 
expressed in a special form below (where c is the speed of light in the E 
frame): 


im 1 / , vdx\ dx v / , vdx\ ,,, dy Irv dz 

dT - ~ y 11 + -J ’ dX = T + a\ df+ ,dY ~ ~d’ dZ ~ ~d’ $$ 

Here we take a Taylor expansion for a, b and d of the form: a~l + om 2 / c 2 + 
0(c~ 4 ); b sw 1 + (3v 2 /c 2 + 0(c~ 4 ); d s=s 1 + Sv 2 /c 2 + 0(c~ 4 ). Recalling that 
7 = ^ 1 2 y 2 from Special Relativity (SR), we see that SR predicts a = — 1/2 
and $ = 1/2. Since SR predicts no contraction in directions orthogonal to a 
boost, it also predicts that (5 = 0. Thus, the RMS parameterizes a possible 
Lorentz violation by a deviation of the parameters (a, /3, 6) from the SR values 
( — 2> 2> 0)- 

By manipulating equation (1) to form the infinitesimals in the S frame, 
we can separate the equation into a boost term (/3 — a — 1), anisotropy term 
{8 —(3+\) and time dilation parameter a+\. Thus, a complete verification of 
LLI in the RMS framework [2, 3] requires a test of (i) the isotropy of the speed 
of light (measuring Pmm = 8 — /3+ 1), a Michelson-Morley (MM) experiment 
[18], (ii) the boost dependence of the speed of light (measuring P KT = /3 — 
a — 1), a Kennedy-Thorndike (KT) experiment [19] and (iii) the time dilation 
parameter (measuring Pis = a+ 1), an Ives-Stillwell (IS) experiment [20, 21]. 
Rotating experiments may be considered Michelson-Morley experiments and 
only measure Pmm, so in this section we restrict ourselves to these types of 
measurements. 

Assuming only a MM type Lorentz violation, and setting ds 2 = <?dT 2 — 
dX 2 — dY 2 — dZ 2 = 0 in E. and transforming according to Eqn. (1) we find 
the coordinate travel time of a light signal in S becomes; 

dt= ^ (p mm x sin 2 0^ + 0(c~ 4 ) (2) 

where dl = \f dx' 2 + dy 2 + dz 2 and 0 is the angle between the Poynting vector 
and the velocity v of S in X. For a modern MM experiment that measures the 
difference frequency between two resonant cavities, the fractional frequency 
difference may be calculated from (2) in a similar way to [11] to give: 

= [/ ( vl *<««>) 2<i<i “ - / ( v -u(®))V] <3) 

Where I j(qj) is the unit vector in the direction of light propagation (Poynting 
vector) of each resonator (labeled by subscripts a and b), and q :] is the variable 
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of integration around the closed path coordinates of the Poynting vector of 
each resonator. 

To calculate the relevant time dependent expressions for v, velocities are 
transformed to a geocentric non-rotating (with respect to distant stars) refer¬ 
ence frame (denoted as the MM-frame) centred at the centre of mass of the 
Earth with its z-axis perpendicular to the equator, pointing north, the x-axis 
pointing towards 11.2h right ascension (aligned with the equatorial projection 
of u defined below). A pictorial representation of the frame is shown in Fig. 1. 
Classical (Galilean) transformations for the velocities and sin 2 # are sufficient 



* 


C 
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Fig. 1. This frame is an Earth centred frame in which the spin axis of the earth is 
the z-axis, and the velocity of the Sun with respect to the CMB is defined to have 
no component in the y direction. Thus, the Earth is spinning at the sidereal rate 
within this frame, and the angle 0 is shown pictorially, but in general is a function 
of position and time as the Earth spins and the experiment rotates. 


as relativistic terms are of order 0(c~ 2 ) and therefore give rise in Eqn. (2) to 
terms of order 0(c~ 5 ). We consider two velocities, the velocity of the sun with 
respect to the CMB u (declination -6.4°, right ascension 11.2 h) and the or¬ 
bital velocity of the Earth v G . Velocities due to the spinning of the Earth and 
laboratory are much smaller and do not impact on the calculations and may 
be ignored. Thus, the sum of the two provide the velocity of the laboratory in 
the universal frame to be inserted in Eqn. (3). In the MM-Earth frame, the 
CMB velocity is: 



u = 


(4) 
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where u « 377km/s and (j)^ ~ —6.4°. To calculate the orbital velocity we first 
consider the Earth in a barycentric non-rotating frame (BRS) with the 2 -axis 
perpendicular to the Earth’s orbital plane and the 2 -axis pointing towards 0° 
right ascension (pointing from the Sun to the Earth at the moment of the 
autumn equinox). 


( -sinAo \ 

cosAo I (5) 

0 J 

where v a « 29.78km/s is the orbital speed of the Earth, and Ao = f2®(t — t Q ) 
with 12® rj 2.0 10 _7 rad/s the angular orbital velocity and t — to the time since 
the autumnal equinox. We first transform to a geocentric frame (GRS) that 
has its x-axis aligned with the BRS one but its 2 -axis perpendicular to the 
equatorial plane of the Earth 


( 6 ) 

where e « 23.27° is the angle between the equatorial and orbital planes of the 
Earth. We then transform to the MM-Earth frame: 

/ v x ^ RS cosa,j, + u^^sino^ \ 
v G = —v x ^ RS sina fJj + t+^coso^ (7) 

V ^ Z o RS ) 

where a^ « 167.9° is the right ascension of u. Summing the two velocities 
from Eqns. (4) and (7) we obtain the velocity of the lab with respect to the 
” universe rest frame”, transformed to the MM-Earth frame 

( ucos</i> M + v 0 (— sinA 0 cosa M + cosAosina^cose) \ 

^(sinAosincc^ + cosAocosa^cose) I. (8) 

iismipn + t' 0 cosAo sine J 

Substituting in the numeric values gives an orbital velocity of (in m/s)-, 

/ 374651 + 5735 cos(A 0 ) + 29118 sin(A 0 ) \ 
v = | —26750 cos(Ao) + 6242 sin(Ao) I . (9) 

\ -42024+ 11765 cos(Ao) J 

The last calculation to make is the time dependence of the the unit vector 
1 along the direction of light propagation, which will depend on the config¬ 
uration of the experiment, including the type of resonator and whether it is 
rotating in the laboratory or not. In section 5 we calculate this dependence 
for a specific experiment, which uses WG modes rotating in the laboratory. 
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2.2 Standard Model Extension 

The SME [6] conglomerates all possible Lorentz-Violating terms and incor¬ 
porates them in a framework, which is an extension of the Standard Model 
of Particle Physics. There are numerous Lorentz-violating terms per parti¬ 
cle sector (i.e. fermions, bosons and photons). However in this work we are 
restricted to the so called minimal “photon-sector”, which only includes 19 
terms. The SME adds additional terms to the Lagrangian of the Standard 
Model for photons. Where as the standard Lagrangian was simply: 

c = (io) 

Under the SME, it becomes [7]: 

C = - \{k F ) KXtlv F KX F^ + \{k AF ) K e KXtlv A x F^ (11) 

where A x is the 4-potential. The (kAF) K terms have the dimensions of mass, 
and are the CPT odd terms [22]. It is argued in [6] that these should be zero 
because they induce instabilities as they are non-negative in the Lagrangian. 
There are also astronomical measurements [7] which place stringent limits on 
kA F ■ From here on these terms are set to zero. 

On the other hand, the ( k F ) Kx „ v terms are CPT even, dimensionless and 
have 19 independent terms out of the 256 possible combinations of k, A, ft. 
and v. Out of these independent Lorentz violating terms, 10 combinations 
have been analysed using astrophysical polarisation tests and have an upper- 
limit of 2 x 10 -32 [7]. This limit is many orders of magnitude less than what 
is expected from laboratory experiments, so these terms are set to zero to 
simplify the calculations and to remain consistent with previous results. 

We can derive the equations of motion for this system by minimising 
the action given by (11), using variational techniques and the definition 
F 1 ' 1 ' = d lt A v — dvAj, and A lt = (<p, A). These equations are similar to those 
of a Maxwellian model in anisotropic media instead of a vacuum. In order to 
express these in a convenient form, we form linear combinations of the CPT 
even term. These are given below [7]: 


= -2{k F ) 0j0k -,K HB ^ = i e^e kra {k F ) pqrs -, 

*DB jk = -*HE kj = {k F ? ivq <M. ( 12 ) 

The dynamics of the model can be described in terms of equivalent B, E, 
H and D fields [7, 11] in a vacuum using the matrices in (12): 

/D\ _ / e °( 1 + K a# \[% KDB \ /E\ 

H \ \f% KHE /‘o X 0 + Kim) J V B / 


(13) 



Resonator Tests of Lorentz Invariance 


Note that (13) is rank 6, as the k matrices are rank 3 as defined in (12). 
The standard Maxwell equations in a vacuum are recovered if these k matrices 
are set to zero. 

Thus the effect of the SME in the photon-sector can be interpreted as 
introducing medium-like properties to the vacuum. In the full SME, this is 
considered as an effect from Planck-scale physics in the early universe. The 
n matrices are all position dependent and thus act as “values” positioned 
throughout space. If one or more of these values are zero, it does not imply 
the rest are also zero as there is no relation between each of the independent 
components. However, there is a linear combination of these components which 
allows us to separate them into birefringent [7] and non-birefringent terms. By 
eliminating those values which have been constrained beyond what we hope 
to achieve in this experiment, these terms can be simply written as in (14). 


(K e+ y k = -KDE ik + (fte—) jfe = ^DE^ ~ KHB^ ~ 


(Ko + ) jk = ~KnB jk + KHE jk ; K- 

Ktr = \s jk {K DE f 



(14) 





Fig. 2. The Sun-Centred Celestial Equatorial Frame (SCCEF), with the Earth 
at different equinoxes during the year. Note that during the vernal equinox, the 
longitude which is at noon has it’s Y axis pointing east, while the longitude which 
is at midnight has it’s Y axis pointing west, and vice-versa for the autumn equinox. 


In the above definitions, K e +, w e _ and are parity-even matrices, while 
k 0 + and k 0 - are the parity-odd matrices. As mentioned in [7, 23], the K e + 
and k 0 - are constrained such that \k jk \ < 2 x 10 -32 . Thus, both n e+ and 
k 0 - are set to zero each time they appear in our equations. Also, k tr is a 
scalar, which resonator experiments are usually insensitive to [24], and is not 
considered in this work. 
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The standard reference frame that we use is the Sun-Centred Celestial 
Equatorial Frame (SCCEF), which is shown in Fig. 2. This is the frame in 
which the sun is at the center, and is inertial with respect to the CMB to first 
order. The axes in this frame are labeled X, Y and Z. The Z axis is defined 
[7, pg. 6] [25, pg. 3] to be parallel to the Earth’s north pole, or 90° declination. 
The X axis points from the sun towards the Earth at the moment of the 
autumn equinox, or 0° right ascension (RA) and 0° declination, while the Y 
axis is at 90° RA and also at 0° declination, usually taken in the J2000.0 
frame. 

The convention described in [7, pg. 18], which has the raised capital indicies 
(J, K) in the SCCEF, has been used. Local coordinates x, y and 2 are defined 
on the Earth’s surface (at the point of the experiment). The 2 axis is defined 
as being locally normal to the ground, vertically upwards. The x axis points 
south and the y axis points east. These coordinates are denoted by the lowered 
capital indicies (j, k) and they rotate with sidereal period AT© = ^ ~ 23 h 
56 min. There is a relation between these two coordinates which is given by 
the following rotation matrix: 

( cosxcosw 0 T e cosxsinw©T 0 -sinx\ 

— sin w 0 T 0 cos w 0 T 0 0 (15) 

sin x cos w 0 T 0 sinxsinu; 0 T 0 sinx / 

Here x is the co-latitude of the laboratory from the north pole, and the T 0 is 
the time coordinate that is related to the sidereal frequency of the Earth. The 
time T 0 is defined in [7] as any time when the y axis and the Y axis align. 
This has been taken to be the first time this occurs after the vernal equinox, 
which points along the negative x-axis. 

When searching for leading order violations it is only necessary to consider 
the spinning of the Earth about itself. However, the orbit of the Earth about 
the sun may also be considered, since it induces Lorentz boosts and we may 
calculate proportional terms to these. Since the Earth moves relatively slowly 
around the Sun compared to the speed of light, the boost terms will be sup¬ 
pressed by the velocity with respect to the speed of light (/3 0 = « 10 -4 ). 

The boost velocity of a point on the Earth’s surface is given by the following 
relation: 

/ sinJ2 0 T \ / —sino; 0 T\ 

/3 = /3© — cost? cos J? 0 T I +/?l cosu> 0 T I (16) 

\ — sin ry cos f2®T J \ 0 / 

Here /3 0 is the value for the boost speed of the orbital motion of the Earth 
and /3l is the boost speed of the lab at the surface of the Earth due to its 
spin motion. The latter is location dependent, but is less than 1.5 x 10 -6 and 
is zero at the poles (y is as defined in Fig. 2). The Lorentz matrix, A(] , that 
implements the transformation from the SCCEF to the laboratory frame with 
the sidereal rotation W J and a boost f3 is given by, 
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/ 1 —/?! -P 2 -P 3 

-(R-P ) 1 R 11 R 12 R 13 

- \-(R-p) 2 R 21 R 22 R 23 
\-{R-p) 3 R 31 R 32 R 33 

After some calculation [7] the k matrices from the SCCEF (indexed by J and 
K) can be express in terms of the values in the laboratory frame (indexed by 
j and k). 

(K DE )i* b = Tt JK {K DE ) JK - T? jk {k db ) jk - T[ kJK {K DB ) JK 
= Tt JK ( KDE ) JK - T^ KJ (K DB ) JK - T( kKJ ( KDB ) JK 
lK DB )t = Tt JK ( KDB ) JK +Tp JK ( KDB ) JK +T{ k - jk ( Khb ) jk 
Where: 

Tl kJK = R jJ R kK 

rpjkJK _ fpPRkJ e KPQpQ 

Here e is the standard anti-symmetric tensor. 

In the section 3 we apply the above to calculate the sensitivity of typical 
resonator experiment. To do this the sensitivity to the components given in 
Eqn. (12) are derived. 


(18) 

(19) 

( 20 ) 


( 21 ) 

( 22 ) 



3 Applying the SME to Resonator Experiments 

The modified Lagrangian of the SME introduces perturbations of the electric 
and magnetic fields in a vacuum. The unperturbed fields are denoted by a 
zero subscript to distinguish them from the Lorentz-violating fields. Putative 
Lorentz violations are produced by motion with respect to a preferred frame, 
which perturbs the fields generating an observable signal. The general frame¬ 
work [7] for denoting the sensitivity of this observable signal in the laboratory 
frame is a linear expression as follows: 

SO = {M DE )iab(KD E yiab + (AIhb^^khb)^ + (MdbJi^kdb)^ (23) 

The summation over the indicies is implied, and the components of the A if k b 
matricies are in general a function of time. The observable is dependent on 
the type of experiment, and in the case of a resonant cavity experiments it 
is the resonance frequency deviation, ^. Since the laboratory frame and the 
resonator frame do not necessarily coincide, we first consider M.{ k s coefficients 
in the resonator frame and later relate it to the laboratory and sun-centred 
frame. 

In general, resonators may be constructed from dielectric and magnetic 
materials. To calculate the M^. k s matricies for such structures a more general 
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form of (13) must be considered, which includes the properties of the medium, 
/.I r (permeability) and e r (permittivity), which are in general second order 
tensors. 

fT>\ / e 0 (e r +KDE) 

V H / ^ Mo 1 (/ x r 1 

Here as was derived in [7], we assume that the fractional frequency shift due 
to Lorentz violations is given by: 


kdb 
+ khb] 


(24) 


[ d 3 x (en 


Av _ 
vo 


1 

“4(h) : 


Here {U) is the energy stored in the field and is given by the standard elec¬ 
trodynamic integral. 


(U) = ^ 3;C (E 0 • Dq + B 0 • Hq) (26) 

In Maxwellian electrodynamics the balance of magnetic and electrical energy 
in a resonator is equal, so the following is true: 

(27) 


This reduces ( U) to an effective normalisation factor for either an electric 
or magnetic filling factor. Also, since the k terms from the integral of (25) 
are only time dependent rather than spatially dependent, the k terms can 
be removed from the integral. Thus, the final term in (25) will be zero since 
the electric and magnetic fields are orthogonal in a resonant structure. By 
equating (23) and (25) in the resonator frame the (AIdb)^ coefficients are 
calculated to be zero, eliminating the possibility of making a measurement 
of Ktr [24, 26]. Assuming the resonator permeability and permittivity have no 
off-diagonal coefficients (i.e. non-gyrotropic) such that; 

fe x 0 0 \ fn x 0 0 \ 

e r = eo I 0 e y 0 I /Ltr- = yu.0 ( 0 fly 0 1 (28) 

\ 0 0 e z ) \ 0 0 n z ) 


the only non-zero coefficients may then be calculated to be 


1 Iv d3x \ E o\ 

{MdeYAs = -- 2 J yd 3 a;Eo *. Eo 


Pe.j 

2i7 


(29) 




Resonator Tests of Lorentz Invariance 


11 


(MHBYrL = Hj 


Iv d3x \ H o\ 

2 f y d 3 xH 0 * ■ H 0 


Pm,j 


(30) 


Thus the M.de and M-hb matricies are diagonal and simply related to the 
electric and magnetic energy filling factors, Pej and Prrij respectively [11]. In 
general a resonator may consist of more than one material, and may include 
vacuum. In this case (29) and (30) may be written more generally (s is the 
number of different materials including vacuum). 


(MdeMs 

ST' Pe j 

(31) 

C MhbYL 

\ ^ /r) Pm) 

~ ^ 2 

(32) 


To measure the resonant frequency it is necessary excite electromagnetic 
fields inside the resonator and then compare it against a similar frequency. To 
be sensitive to violations of LLI, the comparison frequency must be generated 
by a source which exhibits a different dependence on Lorentz violations in 
the photon sector. For example, an atomic standard (such as a hydrogen 
maser) may operate in a mode which is not sensitive to Lorentz violations [11, 
13]. Alternatively, the resonant frequency may be compared against another 
resonator designed to have a different dependence. The latter can be achieved 
by orientating two identical resonators orthogonally [10], or by exciting two 
modes in a matter filled resonator with orthogonal polarizations. In both cases 
the field components must be considered with respect to the laboratory frame 
and not the resonator. For such an experiment the observable becomes the 
frequency difference (between a resonator labeled a and b) such that; 


(33) 


Thus, with respect to the laboratory frame, the effective ( M.DE)iab and 
(. M.HB)iab matricies consistent with (23) become: 


(MDE)a-b = 

/ ( MdE) a X — {-M-De)1 X 0 

0 {Mde)I v - ( Mde)I v 

\ o 0 


0 

0 

(Mde)£ z - ( Mde)I z 


(34) 


(M h 


C Mhb)Z x - ( Mhb)I x 


(35) 


(Mhb)% z ~ ( Mhb)1 


{Mhb)% v - (Mhb)I 
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These equations are general for any resonator experiments, including Fabry- 
Perot and microwave cavity experiments, and simplify the analysis for complex 
resonator configurations, such as whispering gallery mode resonators. Only 
the electric and magnetic filling factors need to be calculated to determine 
the sensitivity coefficients to the observable, which is possible using standard 
numerical techniques [27]. 

To determine the sensitivity of stationary laboratory experiments one cal¬ 
culates the time dependence of (33) due to the sidereal and orbital motion 
of the Earth around the Sun in terms of the Sun-centred coefficients given in 
(14) and (18). This calculation has already been done in [7, 11, 28] and will 
not be repeated here. In the following subsection we generalize this analysis 
to rotating experiments. 


3.1 Rotation in the Laboratory Frame 

Non-rotating experiments [12, 10, 13] that rely on Earth rotation to modulate 
a Lorentz violating effect are not optimal for two reasons. Firstly, the sensi¬ 
tivity is proportional to the noise in the system at the modulation frequency, 
typically best for microwave resonator-oscillators and Fabry-Perot stabilized 
lasers for periods between 10 to 100 seconds. Secondly, the sensitivity is pro¬ 
portional to the square root of the number of periods of the modulation signal, 
therefore taking a relatively long time to acquire sufficient data. Thus, by ro¬ 
tating the experiment the data integration rate is increased and the relevant 
signals are translated to the optimal operating regime [14]. For rotation in the 
laboratory frame the (Al)^ coefficients become a function of time and depend 
on the axis of rotation. In the laboratory it is most practical to rotate around 
the axis of the gravitational field to reduce gravity induced perturbation of the 
experiment. Thus, our analysis includes rotation about the laboratory 2 -axis. 
If we set the time, t = 0, to be defined when the experiment and laboratory 
axes are aligned, and we only consider the time varying components (i.e. the 
most sensitive ones induced by rotation), then for clock-wise rotation of ujr 
rads/sec, (34) and (35)become: 

/ S DE cos(2u R t) -S DE sm(2u R t) 0\ 

{M DE )iab = ~S DE sm(2u R t) -S DE cos(2uj R t) 0 (36) 

\ 0 0 0 / 

Sde = ^ ((M.DE)a X - (Mde)I v - ( Mde)I x + ( Mde)I v ) (37) 

/ S H BCOs(2w R t) -S HB sm(2u R t) 0\ 

(MHB)iab = -Shb sin(2w fi f) -S H b cos(2 u> R t) 0 1 (38) 

\ 0 0 0 / 

Shb = ^ UMhb)T ~ {Mhb)I v ~ ( Mhb)I x + ( Mhb)I v ) 


(39) 
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Note that if one resonator is tested with respect to a stationary generated 
frequency, then the ( M)j J coefficients in the definition of S H b and S H b per¬ 
taining to that frequency must be set to zero. 

To determine the time dependence of the observable (33) we follow the 
same procedure as presented in subsection 2.2 to transform the k matricies 
given in (14) to the k matricies in the laboratory given by (18), (19) and (20). 
We then substitute (37) to (39) into (23) to calculate the time dependence of 
the observable of (33). This is a tedious process and the details are omitted. 
Essentially, because the ( M.DE)iab and (At Hn)iab matricies are time depen¬ 
dent at 2 ujr, the observable signals are at frequencies close to this value and 
are summarized in Table 1. Here the frequency of Earth rotation is defined as 
cu®, and orbit around the sun as 12®. The w® is commonly referred to as the 
sidereal frequency, while the 12® is referred to as the annual frequency. We 
also define the sensitivity factor, S of the experiment as: 

S = Shb — Sde (40) 


u>i 

Cosine Coefficient C Ui /S 

Sine Coefficient S W JS 

2uir 

tsin%)^ 

2/3l sin(x)/5?+ y 

2 ojr + 12© 

sin 2 (x)x 

(2sin(r ? )Kf + y + cos (t/)k* + z ) 

-i/3©sin^(xK + " 

2c or — 12© 

-|/3® sin a (x)x 

(2 sin (r?)/v+ y + cos {r])Ko+) 

|/3© sin 2 {x)m+ 

2u>r + w© 

-2sin^(|)x 

(-PlkZ?+ sin(x)Kt-) 

—2sin^(|)(sin(x)«eA — 0lkI+) 

2 ojr + a;© + 12© 

-2/3© cos(f) sin(r?) sin a (f )k^ 

4/3© cos(|) sin( §) sin 3 (f) x 
(sin(2)^/ +cob(2)«£*) 

2uir + w© — 12© 

-2/3©cos(f)sin(^)sin ;j (fK- 

-4/3©cos(|)cos(|)sin 3 (f)x 

(cos(2) K f/ - sin(2)Kf + z ) 

2 ojr + 2w® 

-sin 4 (|)(^ -Till) 

—2sin 4 (f )kf- 

2ur + 2u><s + 12© 


2/3©sin^(2)sin 4 (f)^ 

2u>r + 2 w© — 12© 

-2/3®cos^(2)sin 4 (f)V + " 

-2/3©cos i (2)sin^)'^ 

2ur — w© 

2cos^(#)(/3^+sin( X )^) 

-2cos*(* )(/3 L m+ +sin(x)KeA) 

2 ujr — a>© + 12© 

wpippir 

-4/3©cos(2)cos 3 (f)sin(|)x 

(cos(2)«f/ -sin(2)«f/) 

2u>r — w© — 12© 

2/3® cos 3 ($) sin(r?) sin( f )kl+ 

4/3©sin(2)cos 3 (|)sin(|)x 

(sin(2)«J/ +C os(2)k^) 

2ujr — 2o>© 

-cos-«") 

2cos 4 (f )ft^_ y 

2ujr — 2o>© + 12© 

- 2/3© cos 2 (2) cos 4 ( § ) K„+ 

2/3© cos^ (^) cos 4 (f) «o+ 

2 lvr — 2a>© — 12© 

2/3© sin y (2) cos 4 (£ )k£+ 

—2/3© sin 2 (2) cos 4 (%) 


Table 1 , Normalized sensitivities with respect to the experiment sensitivity factor 
S for all predicted frequency modulated components. 
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To decorrelate all side bands, more than one year of data is necessary. In 
this case we have eight unknown k coefficients and thirty possible individual 
measurements listed in Table 1, which is an over parameterization. For short 
data sets (less than a year) we do not have enough information to satisfy the 
Nyquist condition to distinguish between frequencies that differ by the annual 
offset (collected in the same blocks). Thus, to make a short data set approx¬ 
imation, we collect the sidebands together (see Fig. 3). The short data set 


Modulation Components vs. Frequency 



Fig. 3. This “frequency stick diagram” shows a schematic of the frequency modu¬ 
lation components of the beat frequency in a convenient form. The sidebands offset 
by ±12® have been trimmed for brevity. Each frequency has two degrees of freedom, 
i.e. a sine and a cosine term, with the phase set by a time t = 0 set with respect to 
the SCCEF. 


approximation is achieved by knowing the angle of the orbit, = l7®t, in the 
sun-centered frame with respect to the negative X-axis (which occurs at the 
vernal equinox as shown in Fig. 2), and then taking a Taylor series expansion 
around that angle. Here we define the phase of the combined rotational and 
sidereal term as 9 and <l>o as the value of <I> when a short data set is taken. 
Since 5$ = $ — o is small with respect to 2ir, via the double angle rule we 
can derive the following relationships: 

sin(0 ± [54> + «P 0 )) = sin(0 ± 5$) cos(^o) ± cos(0 ± 5$) sin(<£ 0 ) 
ss sin(0) cos(??o) ± cos(0) sin(<£ 0 ) 
cos(0 ± (i«P + «P 0 )) = COS {6 ± 5$) cos(^o) T sin(0 ± 5$) sin(£ 0 ) 

*a cos(0) cos(^o) =F sin(0) sin(# 0 ) 


Now we can combine the sidebands as shown in Fig. 3 by applying the above 
relationships to eliminate the dependence on l?®. In this case the components 
from Table 1 decompose to those listed in Table 2. The first feature to notice 
is that the 2a ;r ± 2cjq sidebands are redundant. One might also expect the 
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Coefficient (normalized with respect to S) | 


— cot (^) x S 2 uj r + 2 oj® 


cot*(*) X C 2wH+2wffl 


-2cos^(f)x 

(/3lKo+ + sin (x)(/3® cos(^o ) (cos (r?) - sm(r))K„+) + He-)} 


2cos^(f)x 

(&«?+ + sin( x)(-Ab sin(<A)«? + y + Ab cos(^o) sin (rj)k Y J + kf_ z )) 

s 2wr 

2/3l sin(x)/t^® 

cw-, ( 

-i sin 2 (x) (2/3® cos(^ 0 )(2sin(r?)/to + 1 ' + cos(??)k^ + z ) + 2/3® sin(<P 0 X+ ) 
+ 2sin 2 ( X )«ff) 


-2sitC($)x 

{-0lkZ+ + sin(x)(/3® cos(^ 0 )(cos(r7)K^ + y - sin [n)^o+) + *£- Z )) 

cw We 

-2sin^(f)x 

(-0Lk? + z + sin( X )(-/?® sin(^o)^/ + Ab cos(^o) sin(v)Z z + k? z )) 

<S , 2u> r +2u> (b 

-2sin 4 (f )(/3® sin(^o)K^ + Ab cosfa) cos^ofoY + kfj ) 

C 2ur+2u@ 

- sin 4 (f )(2/3® cosfa) cos(^o)^ - 2/3® sin(<? 0 R + * + (££_* - k\1 )) 


Table 2. Normalized sensitivities with respect to the experiment sensitivity fac¬ 
tor S for all predicted frequency modulated components using the short data set 
approximation. 


2ur ± W 0 sidebands to be redundant as well. The only reason they are not 
is because we have taken into account the velocity of the laboratory due to 
the Earth spinning on its axis, (3 L . In fact it turns out that it is not useful to 
keep this term because /3® is two orders of magnitude larger and when one 
applies the data analysis procedures the sensitivities will be degraded if the 
analysis depends on the (3l terms for the uniqueness of the solution. Thus, 
since it makes no practical sense to keep these terms, we set them to zero. For 
this case the coefficients are listed in Table 3. 

For data sets of less than one year, the components in Table 3 may be used 
to set upper limits on the k coefficients in the SME. Since there are only five 
possible independent components, to set limits on eight coefficients we use 
the same technique as adopted by Lipa et. al. [12]. The k Q+ boost coefficients 
are set to zero to calculate limits on the k e - isotropy coefficients and vice 
versa. This technique assumes no correlation between the isotropy and boost 
coefficients. It would be unlikely that a cancelation of Lorentz violating effects 
would occur, as this would necessitate a fortuitous relationship between the 
coefficients of the same order of value as the boost suppression coefficient (i.e. 
orbit velocity, /3®), and consistent with the correct linear combinations as 
presented in Tab 3. 

Another practical point is that the largest systematic effect occurs at 2 u>r. 
Thus, when setting the limits on the three and coefficients 

we only use the data collected at 2ujr ± w® and 2ujr ± 2 u>® frequencies. 
Likewise for the and (kfj* — k%- ) coefficients. These are 

the same coefficients that have had limits set by the non-rotating experiments 
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Coefficient (normalized with respect to S) 

S2w R -2<jj (r< 

— cot (^) x S 2 uj r + 2 oj® 

C2 UR -2 U „, 

coC(f) x C 2 „ H+2wffl 

s ^)< -V 

coC(f) x S 2uR+ „ m 


-coC(f) x C 2WB+Wffi 

S 2ur 

— 

c 2 „ r 

-■i sin 2 (x) (2/3® cos(^ 0 )(2sin(r?)K^ + y + cos(t])k„+ ) + 2 /?® sin(# 0 )^+) 
+ 2sin 2 ( X )Kf- Z ) 


-2sin 2 (f )(sin( X )(d® cos^oXcosfa)^ - sinfa)*#) + «£-*)) 

C 2UR+Uffi 

-2sin 2 (f)(sin(x)(-/3© S in(^o)^ + 1 ' + 0® cos(^o) MvWo-? +*?-*)) 

&WS+2UJ, 

-2 sin 4 (f )(/3® sin(*o)'# + z + /3® 005(77) cos(# 0 )k^ + ) 

C 2wB + 

- sin 4 (f )(2/3® cos(r?) cos($o)k£+ - 2/3® sin(<5 0 )Ko+ + {*t- ~ )) 


Table 3. Normalized sensitivities with respect to the experiment sensitivity fac¬ 
tor S for all predicted frequency modulated components, using the short data set 
approximation and neglecting components of order (3l- 


[10] [12] [13]. The remaining coefficient *J5(= + %X-) can only be set 

amongst a systematic signal at 2 wr, which is in general much greater than the 
statistical uncertainties at the other frequencies. In this case we can assume all 
coefficients are zero except for the coefficient. However, it is not straight 
forward to set a limit on any putative Lorentz violation amongst a large 
systematic as one can not be sure if the systematic signal actually cancels an 
effect. Since the signal at 2 wr is dominated by systematic effects, it is likely 
that its phase and amplitude will vary across different data sets. In this case 
the systematic signal from multiple data sets can be treated statistically to 
place an upper limit on «f_5. In our experiment we use this technique to set 
an upper of 2.1(5.7) x 10 -14 [15] (see section 6). 

3.2 Phase with respect to the SCCEF 

To extract the k components of the SME out of our observed signal we first 
need to determine the relevant and S$ Wi coefficients listed in Table 1. 
This in turn requires us to know the phase of the experiment’s orientation 
with respect to the SCCEF. In this section we will derive an expression for 
this phase in terms of the time origins of the experiment’s rotation, the Earth’s 
sidereal rotation, and the orbit of the Earth around the Sun. 

In general, we are interested in the frequency components 

2 W[ a ,b] — 2 u)r + aw© + bf2© (41) 

where a and b take on values in the domains 


a£ [-2,2],be [-1,1] 


(42) 
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Thus to determine the b] coefficient we fit the data with a model of the 
form 

cos(2 u> r T r + auj®T® + bf2®T) (43) 

where T)j is the experiment’s rotation time, T 0 is the sidereal time, and T is 
the time since the vernal equinox. 

To simplify our analysis we aim to transform this expression to the form 
cos (at + <j>). To achieve this we note that the difference Sr between the exper¬ 
iment’s rotation time Tr and the time since the vernal equinox T is constant 
over the course of the measurement, as determined by the initial configuration 
of the experiment, and we may write, 

6r = T r - T. (44) 

Similarly the sidereal time and the time since the vernal equinox are related 
by <5®, 

<5 0 = T 0 - T (45) 

By combining (43), (44) and (45) we arrive at an expression of the desired 
form. 


cos(2 ujrTr + au>®T® + bf2®T) = cos(2wfj(cif{ + T) + aw 0 (5 0 + T) + bf2®T) 
= cos((2w fl + aw 0 + bfi®)T 
+2 ujrSr + aw®6®) 

Thus we can account for the phase of the experiment relative to the SCCEF 
by determining Sr and <5 0 . The origin of the experiment’s rotation time Tr 
is defined to be the instant at which the axis of symmetry of the first res¬ 
onator (resonator a) is aligned with the local y axis. Our experiment has been 
designed such that the time origin of the data acquisition coincides with the 
same event, rendering Sr = 0 in our case. 

We also need to obtain <J 0 for the sidereal rotation. We define T 0 = 0 as 
in [7] to be the instant the local y axis and the SCCEF Y axis are aligned 
(noon) in the laboratory (see figure 2). Let us define T v to be the time in 
seconds after midnight UTC+0, at which the vernal equinox has occurred 
in the J2000.0 frame [7]. For convenience we also define our longitude Tl 
in terms of sidereal seconds from midnight (in the case of our laboratory 
Ti = 115.826° x 23h 3QQ™ m = 27721sec). There exists a special location whose 
meridian is at noon at the vernal equinox. For this special location (during 
the vernal equinox), £ 0 = 0 since the time when the y and Y axes align and 
the vernal equinox are the same. We see geometrically that any longitude 
greater than this meridian will have positive 6®, otherwise if the longitude is 
less than this meridian it would have negative 6®. As shown in Fig. 4, we can 
now derive an expression for <5 0 . 


6® =T l + T v 


23hr56min 


(46) 


Hence we are able to determine the phase of the experiment’s orientation 
relative to the SCCEF. 
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Y 

4 


Midnight 



UTC+0 


Noon 


Fig. 4. Diagram showing meridians and angles. 


4 Comparison of sensitivity of various resonator 
experiments in the SME 

In this section we show how the general analysis may be applied to some 
common resonator configurations for testing LLI. Also, we propose some new 
configurations based on exciting two modes in matter filled resonators. The 
comparison is made by calculating the sensitivity parameter S of the resonator 
using Eqns. (31) to (40). Note that the sign of the S factor depends on the 
definition of the first resonator. Practically this will need to be the resonator 
that exhibits the largest value of frequency. In this work, where appropriate, 
we assume the first resonator is aligned along the y-axis. 

4.1 Fabry-Perot resonators 

Experiments based on laser stabilized Fabry-Perot resonators typically use 
either one [8] or two [10] cavities placed with the lengths orthogonal to the 
laboratory 2 -axis. In a vacuum filled cavity it is easy to show that \S\ = | 
for the configuration in Fig .5. In contrast, when one rotating cavity is com¬ 
pared to a stationary one the value is reduced by a factor of 2, to |«S| = \. 
It is also interesting to consider the sensitivity of matter filled cavities in the 
photon sector. Here, for simplicity we assume the relative permeability and 
permittivity are scalars of jj r and e r respectively. It is straight forward to 
add anisotropy and only modifies the sensitivity slightly, so for brevity is not 
considered here. If similar configurations to Fig. 5 are constructed from solid 
material the sensitivity factor, S, becomes dependent on polarization. This 
effect also allows for a sensitive experiment by exciting two modes of different 
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Fig. 5. Typical configuration of a rotating Lorentz invariance test using Fabry-Perot 
cavities. 


polarization inside one cavity (Dual-Mode), of which some examples are shown 
in Fig . 6. Such cavities have been built previously to measure birefringence 
[29]. High finesse matter cavities can be made by using low-loss crystalline 
dielectric materials at optical frequencies [30, 31]. The sensitivity for different 
configurations are compared in Table 4. For a low-loss dielectric material with 



Fig. 6. New proposed matter filled Fabry-Perot cavity configurations in which two 
modes of orthogonal polarization are excited and compared, and are sensitive to 
violations in Lorentz invariance in the photon sector of the SME 


E z polarization in the two orthogonal cavities (Fig .5) the sensitive factor, S, 
is the same as the vacuum cavity, while for the circularly polarized case, the 
sensitivity is close to that of the single vacuum cavity resonator. In contrast 
the same experiment with H z polarization has reduced sensitivity of the or- 
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Configuration 

Sensitivity Factor S 

Fig . 5 E z 


Fig . 5 H z 


Fig . 5 Circular Polarization 

i (£+*■) 

Fig . 6 (a) 


Fig . 6 (b) 



Table 4. Value of the S factor for various configurations of Fabry-Perot cavity 
experiments 


der of the permittivity of the material. The sensitivity of the two Dual-Mode 
resonators gives the possibility of realizing a similar sensitivity to dual cav¬ 
ity experiments, but within the same cavity. The configuration should have 
a large degree of common mode rejection, and will be much more insensitive 
to external effects like temperature, vibration etc. and other systematics, and 
may be worth pursuing for these reasons. Note that Muller has recently com¬ 
pleted an analysis of conventional cavity configurations in the electron (due 
to dispersion changes) and photon sector [32]. In our analysis we have only 
considered the photon sector and we have proposed some new unconventional 
configurations. It may be interesting to analyze these configurations in the 
electron sector. In the next subsection we consider similar configurations for 
Whispering Gallery (WG) modes. 


4.2 Whispering gallery mode resonators 

In this subsection we consider ’pure’ WG modes, with the electric and mag¬ 
netic fields propagating around with cylindrical symmetry at a discontinuity, 
with the direction of the Poynting vector (E x B) as shown in figure 7. Thus, 
it is natural to analyse such modes in cylindrical coordinates {r, <f), zj. 

For an actual WG mode the wave is reflected off an electromagnetic dis¬ 
continuity, and the fields mainly lie within the radius of the discontinuity and 
a smaller inner caustic[ll]. However, by taking the limit as the azimuthal 
mode number m tends to infinity, the inner caustic converges to the radius of 
the discontinuity and the fields are reduced to a Dirac delta function. There 
are two possible polarizations, WGE with dominant H z and E r fields and 
WGH with dominant E z and H ,. For ’pure’ WG modes, WGE have non-zero 
electric and magnetic filling factors of Pe r = 1 and Pm z = 1, and WGH have 
electric and magnetic filling factors of Pe z = 1 and Pm r = 1. in cylindrical 
coordinates. The electric and magnetic filling factors may be converted from 
cylindrical to cartesian symmetry by (the z component of the filling factor 
need not be transformed): 


Pe x = Pe l 


Pe r + Pe q 


Pm x = Prrij 


Pm r + Pm q 


(47) 
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We can now do a similar analysis to subsection 4.1 for configurations shown 
in Fig. 8 and 9 for the WG case, with the computed sensitivities listed in Table 
5. In vacuum the S factor is half that of the FP cavities in subsection 4.1, 
and the Dual-Mode resonator is insensitive. However, in a low loss dielectric 
the S factor approaches the same value for WGE modes as the FP cavity 
experiments, but the WGH modes remain about a factor of two less sensitive. 
The value of the S factor for the Dual-Mode resonator is the mean value of 
the WGE and WGH modes. 
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Fig. 9. Rotating Lorentz invariance test using two WG modes excited in one cavity 
(Dual-Mode resonator) 



Table 5. Value of the S factor for various configurations of WG mode resonator 
cavity experiments 


We have shown that similar sensitivities can be achieved with FP and WG 
cavity resonators. At UWA we have developed an experiment that uses low 
loss sapphire crystals, which exhibit a small uniaxial dielectric anistropy. The 
calculations of the sensitivity are presented section 6. 

5 Applying the RMS to Whispering Gallery Mode 
Resonator Experiments 

In this section we restrict ourselves to analysis of whispering gallery mode 
resonator experiments, as the analysis has been well described for Fabry- 
Perot resonators previously [28]. For the whispering gallery mode experiment 
as shown in Fig. 8, the variable of integration around the path of the resonator 
is naturally chosen as the azimuthal angle, 4>j, relative to the cylindrical co¬ 
ordinates of each resonator. Thus, from Eqn. (3) a frequency shift due to a 
putative Lorentz violation in the RMS framework is given by, 

y? = [/ - / (v-Mw) 2 


(48) 




Resonator Tests of Lorentz Invariance 


23 


The dominant components of the velocity vector v were already calculated 
in section 2.1, so to complete the calculation the time dependence of I„ and 
If, must be calculated with respect to the MM-Earth frame. This of course 
depends on the sidereal and semi-sidereal frequencies, as well as the rotation 
frequency of the experiment. To start the calculation we define the time, 
t = 0 when the axis of the two WG resonators are aligned as shown in Fig. 
8 (i.e. the resonators align with the laboratory frame). Then from this time 
we assume the resonator is rotated in a anti-clockwise direction of frequency 
uj a , so the angle of rotation is 7 = u s (t — t s ). Also, the longitudinal angle of 
the experiment is A, which is dependent on the sidereal frequency and given 
by A = w®(i — ti). Then we define the resonator with its cylinder axis in the 
y direction as resonator a, and the resonator with its cylinder axis in the x 
direction as resonator b. We also assume the WG modes are oscillating in 
a clockwise direction. In actual fact the calculation has been verified to be 
independent of the WG mode direction, and in most experiments is usually 
a standing wave (depending on the excitation) [33]. Thus in the laboratory 
frame at t = 0 the unit vectors in the direction of the Poynting vector are; 



/ — sin <j> a \ 

latta) = 0 

( ° \ 

I b(<t>b) = COS(/) b 1 

(49) 


\ COS (j> a ) 

V sin</>& / 


Now if we 

transform from the resonator to the laboratory, then to the 

MM- 

Earth frar 

ne the unit vectors become. 




/ —sin0 a (cosAcosxcosq 

— surysinA) + cos0 a cosAsinx\ 


I Eartl 

\:a= —sin^cosAsnry + co 

sqcosxsinA) + cos</> a sinAsinx I 

(50) 


\ cosycosd,, 

; + sinxcos7sin</> a J 



/ — sin<^(,(cosAcosxsin7 

+ cosqsinA) + cos</>&cosAsinx \ 


I Ear 

th-.b = sin</>;,(cosAcos7 - sii 

iqcosxsinA) + cos^sinAsinx 1 

(51) 


\ cosxcos c/)b + sinxsiirysini/)*, J 


Here as ii 

1 the previous sections x is 

the angle from the north pole (co- 


latitude). 

The next step is to substitute (50), (51) and (9) into (48). However, to be 
consistent with the SME analysis the phase should be calculated with respect 
to the vernal equinox, so that Ao = <&o + n is substituted into Eqn. (9) before 
we substitute it into (48) to calculate the frequency shift. Also, because we 
defined the rotation to be clockwise in the SME, to be consistent we define 
7 r = uJn(t — tg) where ujr = — ui s . In this case the frequency components, 
which experience a frequency shift are given in Table 6. From the results of 
the calculation we note that perturbations due to Lorentz violations occur at 
the same frequencies as the SME (see subsection 7.2). Fortunately, it is not 
necessary to consider perturbations at exactly twice the spin frequency, 2u s , 
that are primarily due to the larger systematic effects associated with the 
rotation, as we only need to put a limit on one parameter. Also, the cosine 
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Ui 

10 'CuvJPmm 

2 ojr 

-3.904 + 3.904cos(x) 2 + 0.098sin(x) 2 
— sin(^o) (—0.607 + 0.607cos(x) 2 ) 

- cos(^o) (-0.120 + 0.120cos(x) 2 - 0.055sin(x) 2 ) 

2 ojr + o>© 

-0.876 sin(x) + 0.876 cos(x) sin(x) 

+ sin(<A,) (0.068 sin(x) - 0.068 cos(x) sin( X )) 

+ cos(^o) (-0.232 sin(x) + 0.232 cos(x) sin( X )) 

2u>r — o»® 

0.876 sin(x) + 0.876 cos(x) sin(x) 

+ sin(£o) (-0.068 sin(x) - 0.068 cos(x) sin(x)) 

+ cos(^o) (0.232 sin(x) + 0.232 cos( X ) sin( X )) 

2c jr + 2u>© 

1.952 - 3.904 cos(x) + 1.952cos(x) 2 

Tcos(4*o) (-0.060 + 0.120 cos(x) - 0.060cos(x) 2 ) 
+ sin(<A,) (-0.303 + 0.607cos(x) - 0.303cos(x) 2 ) 

2u)r — 2u>© 

1.952 + 3.904 cos(x) + 1.952cos(x) 2 
+ cos(^ 0 ) (-0.060 - 0.120 cos(x) - 0.060cos(x) 2 ) 
+ sin(<£ 0 ) (-0.303 - 0.607cos(x) - 0.303cos(x) 2 ) 


Table 6. Dominant coefficients in the RMS, using a short data set approximation 
calculated from Eqn. (48). 


components ( Cu Ui ) with respect to the CMB are the most sensitive, so we 
need not consider the sine components. 

5.1 Phase with respect to the CMB 

To extract the Pmm term from our data we must first determine the phase of 
our experiment with respect to the CMB. Thus, in similar way to the reasoning 
for the SME (see subsection 3.2) we require Sr, the difference between the 
experiment’s rotation time and the time since the vernal equinox, and (5®, the 
difference between the sidereal time and the time since the vernal equinox. 

As was the case for the SME, Sr = 0 since the axis of symmetry of the 
first resonator, a, is aligned with the local y-axis at Tr = 0. However, 5® 
will be different since in the case of the RMS it is measured with respect 
to the CMB (or MM-Earth frame), not the SCCEF. The CMB is oriented 
at 11.2 h right ascension, 6.4 degrees declination relative to the equatorial 
plane. Let us define T v to be the time in seconds after midnight UTC+0, 
at which the vernal equinox has occurred in the J2000.0 frame [7]. T u is the 
direction of the CMB (11.2h). For convenience we also define our longitude 
T l in terms of sidereal seconds from midnight (in the case of our laboratory 
Ti = 115.826° x m 36Qo min = 27721sec). As shown in Fig. 10, we now have an 
expression for 5®. 


<5® = T l + T v - (T u + 


23hr56min , 


(52) 
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Fig. 10. Diagram showing meridians and angles used to determine the phase of the 
experiment with respect to the CMB. 


Hence we are able to determine the phase of the experiment orientation 
relative to the CMB. 


6 The University of Western Australia Rotating 
Experiment 

Our experiment consists of two cylindrical sapphire resonators of 3 cm diam¬ 
eter and height supported by spindles at either end within superconducting 
niobium cavities [34], which are oriented with their cylindrical axes orthog¬ 
onal to each other in the horizontal plane (see Fig. 11). Whispering gallery 
modes [35] are excited close to 10 GHz, with a difference frequency of 226 kHz. 
The frequencies are stabilized using a Pound locking scheme, and amplitude 
variations are suppressed using an additional control circuit. A detailed de¬ 
scription of the cryogenic oscillators can be found in [36, 37], and a schematic 
of the experimental setup is shown in Fig. 12. The resonators are mounted in 
a common copper block, which provides common mode rejection of tempera¬ 
ture fluctuations due to high thermal conductivity at cryogenic temperatures. 
The structure is in turn mounted inside two successive stainless steel vacuum 
cylinders from a copper post, which provides the thermal connection between 
the cavities and the liquid helium bath. A stainless steel section within the 
copper post provides thermal filtering of bath temperature fluctuations. A 
foil heater and carbon-glass temperature sensor attached to the copper post 
controls the temperature set point to 6 K with mK stability. Two stages of 
vacuum isolation are used to avoid contamination of the sapphire resonators 
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Fig. 11. The two resonators are positioned orthogonal to each other in the mounting 
structure. One of the sapphires can be seen mounted inside the superconducting 
niobium cavity. The spindles are firmly held in each lid by sprung brass bushes. 


from the microwave and temperature control devices located in the cryogenic 
environment. 

A schematic of the rotation system is shown in Fig. 13. A cryogenic dewar 
containing the resonators, along with the room temperature oscillator circuits 
and control electronics, is suspended within a ring bearing. A multiple ”V” 
shaped suspension made from loops of elastic shock cord avoids high Q-factor 
pendulum modes by ensuring that the cord has to stretch and shrink (pro¬ 
viding damping losses) for horizontal motion as well as vertical. The rotation 
system is driven by a microprocessor controlled stepper motor. A commercial 
18 conductor slip ring connector, with a hollow through bore, transfers power 
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and various signals to and from the rotating experiment. A mercury based 
rotating coaxial connector transmits the difference frequency to a stationary 
frequency counter referenced to an Oscilloquartz oscillator. The data acquisi¬ 
tion system logs the difference frequency as a function of orientation, as well 
as monitoring systematic effects including the temperature of the resonators, 
liquid helium bath level, ambient room temperature, oscillator control signals, 
tilt, and helium return line pressure. 

Inside the sapphire crystals standing waves are set up with the dominant 
electric and magnetic fields in the axial and radial directions respectively, 
corresponding to a propagation (Poynting) vector around the circumference. 
The observable of the experiment is the difference frequency, and to test for 
Lorentz violations the perturbation of the observable with respect to an alter¬ 
native test theory must be derived. The mode which we excite is a Whispering 
Gallery mode we have a choice of WGE m . Ti;P or WGH mi „ iP modes, the first 
subscript, m, gives the azimuthal mode number, while n and p give the num¬ 
ber of zero crossings in the radial and z-direction respectively. Typically the 
so called fundamental mode families WGE mj o,o or WGH TO) o,o as they have the 
highest Q-factors. To calculate the sensitivity in the RMS we use the tech- 
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Fig. 13. Schematic of the cryogenic dewar, mounted in the rotation table. Inside 
the dewar a schematic of the two orthogonally orientated resonators is shown, along 
with the Poynting vectors of propagation Si and S2 ■ 


nique presented in section 5, while in the following subsection we numerically 
compute the sensitivity in the SME. 

6.1 Sensitivity in the SME 

In this subsection we calculate the sensitivity of the fundamental WG mode 
families, WGE mi o,o and WGH TO) o,o to putative Lorentz violation in the SME, 
and compare it with the ’pure’ WG approximation given in Fig. 8. For a 
proper analysis of the sapphire loaded cavity resonators two regions of space 
need to be taken into account: the anisotropic crystal and the cavity free 
space surrounding it (see Fig. 14). The latter has a relative permittivity of 1, 
while both have relative permeability of 1 in all directions. The calculations 
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proceed by splitting up V into V\ (the crystal) and V 2 (the freespace), so we 
may sum the components of the M. matricies over the two volumes (see Eqns. 
(31) and (32)). The resonator operates close to liquid helium temperatures (6 
Kelvin), where the permittivity of sapphire is, ej_ = 9.272 and ey = 11.349. To 



Fig. 14. Schematic of a cylindrical sapphire crystal resonator, with dimensions in 
mm. The crystal exhibits uniaxial anisotropy with the axis of symmetry (c-axis) 
aligned along the cylindrical z axis. The permittivity along the c-axis given by e\ \. 
Perpendicular to the c-axis in the x, y or r, <j> plane, the permittivity is given by ex- 
The two regions shown are 1. the crystal and 2. the cavity. 


determine the sensitivity, we need to just calculate the experiments S factor 
in a similar way to the calculation for the ’pure’ WG modes in subsection 
4.2. In this case the electric and magnetic filling factors must be calculated 
using a numeric technique such as finite element analysis, method of lines or 
separation of variables [11]. In this work we have chosen to use method of lines 
developed at IRCOM at the University of Limoges [39] . The calculated electric 
and magnetic field densities for the chosen mode (WGHs : o,o) of operation at 
10 GHz is shown in figure 15, and the S factor is calculated to be 0.19575. 

The actual WG modes have all field components in both regions of the 
crystal. This modifies the sensitivity slightly, but approaches the limit of the 
’pure’ WG mode as m —> 00 . The magnitude of the S factor for the funda¬ 
mental WGE and WGH modes at X-Band (8GHz-12GHz) are plotted in Fig. 
16. The WGH modes seem converge nicely towards the predicted ’pure’ WGH 
mode sensitivity, while the WGE modes have a dip in sensitivity. This can 
be explained by an intersection with another mode of the same m number, 
resulting in a spurious mode interaction[38] . This does not occur in WGH 
modes since they are the lowest frequency modes for the mode number m 
(refer to figure 2 of [38]). It is important to note that about a factor of two in 
sensitivity can be gained if we use a WGE mode rather than a WGH mode. 
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Electric r t 



Magnetic r <]> 



Fig. 15. The magnetic and electric field density plots shown in the top right-hand 
quadrant of an axial slice through the sapphire crystal and cavity for the WGHs, 0,0 
mode. Note the dominant fields are E z and H, consistent with a pure whispering 
gallery mode approximation. 


However because we are using a 3 cm crystal rather than a 5 cm crystal, the 
Q-factor of WGE modes are degraded due to radiation and wall losses. In 
the future we can markedly improve the sensitivity by employing the typical 
5 cm cavities that operate in WGE modes, as were used in the non-rotating 
experiments [11]. 


7 Data Analysis and Interpretation of Results 

Fig.17 shows typical fractional frequency instability of the 226 kHz differ¬ 
ence with respect to 10 GHz, and compares the instability when rotating and 
stationary. A minimum of 1.6 x 10 -14 is recorded at 40s. Rotation induced 
systematic effects degrade the stability up to 18s due to signals at the rotation 
frequency of 0.0567/2 and its harmonics. We have determined that tilt varia¬ 
tions dominate the systematic by measuring the magnitude of the fractional 
frequency dependence on tilt and the variation in tilt at twice the rotation 
frequency, 2ujr{{).\1Hz), as the experiment rotates. We minimize the effect of 
tilt by manually setting the rotation bearing until our tilt sensor reads a min¬ 
imum at 2u>r. The latter data sets were up to an order of magnitude reduced 
in amplitude as we became more experienced at this process. The remaining 
systematic signal is due to the residual tilt variations, which could be further 
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Fig. 16. The frequency and |<S| factor as a function of mode number, m for the two 
lowest frequency WGE and WGH mode families. Note the general convergence to 
the ’pure’ WG value as m increases. 



Fig. IT. Square Root Allan Variance fractional frequency instability measurement of 
the difference frequency when rotating (crosses) and stationary (circles). The hump 
at short integration times is due to systematic effects associated with the rotation 
of the experiment, with a period of 18 seconds. Above 18 seconds the instability is 
the same as when the experiment is stationary. 


annulled with an automatic tilt control system. It is still possible to be sen¬ 
sitive to Lorentz violations in the presence of these systematics by measuring 
the sidereal, Uq, and semi-sidereal, 2w 0 , sidebands about 2uj r , as was done 
in [8]. The amplitude and phase of a Lorentz violating signal is determined 
by fitting the parameters of Eq. (53) to the data, with the phase of the fit 
adjusted according to the test theory used. 
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- = A + Bt + ^2 C Ui cos(u)it + ipi) + S Ui sm(w$ + ipi) (53) 

Here i/q is the average unperturbed frequency of the two sapphire resonators, 
and Avq is the perturbation of the 226 kHz difference frequency, A and B 
determine the frequency offset and drift, and C Ui and S uli are the amplitudes 
of a cosine and sine at frequency w* respectively. In the final analysis we 
fit 5 frequencies to the data, w* = (2u}r,2ujr ± w ®, 2ujr ± 2<j®), as well as 
the frequency offset and drift. The correlation coefficients between the fitted 
parameters are all between 10 -2 to 10 -5 . Since the residuals exhibit a signif¬ 
icantly non-white behavior, the optimal regression method is weighted least 
squares (WLS) [13]. WLS involves pre-multiplying both the experimental data 
and the model matrix by a whitening matrix determined by the noise type of 
the residuals of an ordinary least squares analysis. 

We have acquired 5 sets of data over a period of 3 months beginning 
December 2004, totaling 18 days. The length of the sets (in days) and size 
of the systematic are (3.6,2.3 x 10 -14 ), (2.4,2.1 x 10 -14 ), (1.9,2.6 x 10 -14 ), 
(4.7,1.4x 10 -15 ), and (6.1,8.8 x 10 -15 ) respectively. We have observed leakage 
of the systematic into the neighboring side bands due to aliasing when the 
data set is not long enough or the systematic is too large. Fig. 18 shows the 
total amplitude resulting from a WLS fit to 2 of the data sets over a range of 
frequencies about 2 u>r. It is evident that the systematic of data set 1 at 2 ivr 
is affecting the fitted amplitude of the sidereal sidebands 2 u>r ± w® due to its 
relatively short length and large systematics. By analyzing all five data sets 
simultaneously using WLS the effective length of the data is increased, reduc¬ 
ing the width of the systematic sufficiently as to not contribute significantly 
to the sidereal and semi-sidereal sidebands. 


7.1 Standard Model Extension Framework 

In the photon sector of the SME 10 independent components of k e+ and k a - 
have been constrained by astronomical measurements to < 2 x 10 -32 [7, 23]. 
Seven components of k e - and k a+ have been constrained in optical and mi¬ 
crowave cavity experiments [10, 13] at the 10 -15 and 10 -11 level respectively, 
while the scalar k tr component recently had an upper limit set of < 10 -4 [24]. 
The remaining /tpp component could not be previously constrained in non¬ 
rotating experiments [10, 13]. In contrast, our rotating experiment is sensitive 
to kfz. However, it appears only at 2 ujr, which is dominated by systematic 
effects. By using the formulas derived in Table 3 for short data sets and the 
S factor for the IFGL/go.o mode in Fig. 16, the resulting numerical relation 
between the parameters of the SME and the C Wi and S Ui coefficients were 
calculated and are given in Table 7. 

From our combined analysis of all data sets, and using the relation to 
given in Table 7, we determine a value for kf^ of 4.1(0.5) x 10 -15 . However, 
since we do not know if the systematic has canceled a Lorentz violating signal 
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Frequency - 2co R 

Fig. 18. Spectrum of amplitudes \JCS, i + calculated using WLS, showing sys¬ 
tematic leakage about 2u>r for 2 data sets, data set 1 (3.6 days, circles), data set 5 
(6.1 days, squares) and the combined data (18 days spanning 3 months, solid trian¬ 
gles). Here w® is the sidereal frequency (11.6 ^Hz). By comparing a variety of data 
sets we have seen that leakage is reduced in longer data sets with lower systemat- 
ics. The insets show the typical amplitude away from the systematic, which have 
statistical uncetainties of order 10 -16 . 


at 2 Ur, we cannot reasonably claim this as an upper limit. Since we have five 
individual data sets, a limit can be set by treating the C 2wh coefficient as a 
statistic. The phase of the systematic depends on the initial experimental con¬ 
ditions, and is random across the data sets. Thus, we have five values of C^r, 
({—4.2,11.4,21.4,1.3, —8.1} in 10 -15 ), two are negative coefficients and three 
are positive. If we take the mean of these coefficients, the systematic signal 
will cancel if the phase is random, but the possible Lorentz violating signal 
will not, since the phase is constant. Thus a limit can be set by taking the 
mean and standard deviation of the five coefficient of C^uir • This gives a more 
conservative bound of 2.1(5.7) x 10 -14 , which includes zero. Our experiment 
is also sensitive to all other seven components of k e - and k a+ (see Table 7) 
and improves present limits by up to a factor of 7, as shown in Table 8. 

7.2 Robertson, Mansouri, Sexl Framework 

From Eqn. (48), the dominant coefficients are calculated to be only due to the 
cosine terms with respect to the CMB right ascension, Cu Ui , and the theory 
predicts no perturbations in the quadrature term. Since our experiment is 
rotating clock wise we can substitute 2 u> s = —2ujr, and once we perform 
the integral and substitute all the numeric values. Following this method we 
calculate the coefficients as shown in Table 9. 

The same five data sets were then re-analysed in the correct quadrature 
with respect to the CMB, with the results listed with the coefficients in Ta¬ 
ble 9. The measured and statistical uncertainty of Pmm is determined to be 
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Wi 

c Ui 

S Mi 

2u>r 

0.21 ftf- 


2 OJR + w® 

2.5 x 10 _6 sin<£ 0 Ko+ 

— 1.0 x 10 -5 cos^ok^h? 
-0.27 k£_ z 

-COS &0 [2.3 x 10~ 5 ko+ - 1.0 x 10“ 5 /to+J 
-0.27 kX- 

2 u)r + 2o>® 

-2.1 x 10“ 5 cos $ 0 kf+ 
+2.3 x 10 -5 sin^oK^n? 
-O.U(k™-kD 

-2.3 x 10 -5 sin <&ok*+ 

—2.1 x 10 -6 cos^ok^h? — 0.23 k*_ 

2u)r — w® 

—0.31C2 WH +u; ffi 

O.31$2a>#+2u>0 

2u>r — 2o>® 

9.4 X 10- a C'2u, Jl +2a, (B 

-9.4 x 10-^ 2wjf+2wm 


Table 7. Coefficients C Ui and S Wi in (1) for the five frequencies of interest and their 
relation to the components of the SME parameters k e - and k 0 +, derived using a 
short data set approximation including terms up to first order in orbital velocity, 
where 4>o is the phase of the orbit since the vernal equinox. Note that for short 
data sets the upper and lower sidereal sidebands are redundant, which reduces the 
number of independent measurements to 5. To lift the redundancy, more than a year 
of data is required so annual offsets may be de-correlated from the twice rotational 
and sidereal sidebands listed. 





ke- 

(*t- ~ He- ) 

this work 

-0.63(0.43) 

0.19(0.37)1 

-0.45(0.37)1 

-1.3(0.9) 

from [13] 

-5.7(2.3) 

-3.2(1.3) 

-0.5(1.3) 

-3.2(4.6) 


Ke- 

Ko V 

k„+ 

K+ 

this work 

21(57) 

0.20(0.21)| 

-0.91(0.46)| 

0.44(0.46) 

from [13] 

- 

-1.8(1.5) 

-1.4(2.3) 

2.7(2.2) 


Table 8. Results for the SME Lorentz violation parameters, assuming no cancelation 
between the isotropy terms. k e - (in 10 -15 ) and first order boost terms k 0 + (in 10 -11 ) 
[12]- 


m 

Cu Ui 

Pmm 

2uir + w® 

[-1.13 x 10 ' - 3.01 x 10~ 8 cos$ 0 
+8.83 x 10- 9 sin<£ 0 \P M m 

—2.1(7.2) 

2 U1R — w® 

[3.51 x 10 -8 + 9.31 x 10" 9 cos 
—2.73 x 10- 9 sin^ o ]P M m 

62.4(23.3) 

2 ojr + 2u>® 

[4.56 x 10“ 7 - 1.39 x 10“ 8 cos^o 
-7.08 x 10 -8 sin 4>o\Pmm 

—1.3(2.1) 

2ujr — 2w® 

[4.37 x 10“ 8 - 1.34 x 10“ 9 cos^o 
-6.78 x 10- 9 sin^ o ]P M m 

-7.5(22.1) 


Table 9. Dominant coefficients in the RMS, using a short data set approximation 
calculated from Eq. (48). The measured values of Pmm (in 10 -10 ) are shown together 
with the statistical uncertainties in the bracket. From this data the measured and 
statistical uncertainty of Pmm is determined to be —0.9(2.0) x 10 -10 , which repre¬ 
sents more than a factor of 7.5 improvement over previous results 2.2(1.5) x 10 -9 [10]. 
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—0.9(2.0) x 10 _1 °, which represents a factor of 7.5 improvement over previous 
results 2.2(1.5) x 1CT 9 [10]. 


8 Summary 

Rotating resonator experiments are emerging as one of the most sensitive types 
of Local Lorentz Invariance tests in electrodynamics (also see other contribu¬ 
tions within these proceedings). In this work we have analysed in detail such 
experiments to putative Lorentz violation in both the RMS and SME frame¬ 
works. In the RMS, rotating experiments only enhance the sensitivity to the 
Michelson-Morley parameter, P MJf , and are not sensitive to the Kennedy- 
Thorndike , Pkt, or Ives-Stilwell, Pis , parameters. In the SME non-rotating 
resonator experiments in the laboratory test for four components of the k e - 
tensor and three components of the k a+ tensor, with the scalar coefficient k tr 
and unmeasurable. Rotation in the SME enhances the sensitivity to the 
seven components and also allows the determination of the component. 
We have shown that all resonator experiment exhibit the same relative fre¬ 
quency spectrum to the putative signal to within a multiplicative sensitivity 
factor, S. This was utilized to compare the sensitivity of different FP and 
WG resonator configurations, leading to the proposal of some new dual-mode 
resonator experiments. 

We applied the above analysis to our experiment at the University of 
Western Australia, which is based on rotating cryogenic sapphire whispering 
gallery mode microwave oscillators. In summary, we presented the first results 
of the experiment, which we set bounds on 7 components of the SME photon 
sector (Table 8) and Pmm (Table 9) of the RMS framework up to a factor of 7.5 
more stringent than those obtained from previous experiments. We also set an 
upper limit (2.1(5.7) x 10 -14 ) on the previously unmeasured SME component 
~zz q-Q f ur ther improve these results, tilt and environmental controls will be 
implemented to reduce systematic effects. To remove the assumption that the 
k a+ and k e - do not cancel each other, data integration will continue for more 
than a year. 
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